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ABSTRACT 

We study clusters in the BCS cluster sample which are observed by Chandra and are more 
distant than redshift, z > 0.1. We select from this subsample the clusters which have both a 
short central cooling time and a central temperature drop, and also those with a central radio 
source. Six of the clusters have clear bubbles near the centre. We calculate the heating by these 
bubbles and express it as the ratio rhcat Acooi = 1.34 ± 0.20. This result is used to calculate 
the average size of bubbles expected in all clusters with central radio sources. In three cases 
the predicted bubble sizes approximately match the observed radio lobe dimensions. 

We combine this cluster sample with the B55 sample studied in earlier work to increase 
the total sample size and redshift range. This extended sample contains 71 clusters in the 
redshift range ^ z ^ 0.4. The average distance out to which the bubbles offset the X- 
ray cooling in the combined sample is at least fheat/fcooi = 0.92 ± 0.11. The distribution 
of central cooling times for the combined sample shows no clusters with clear bubbles and 
icooi > 1-2 Gyr. An investigation of the evolution of cluster parameters within the redshift 
range of the combined samples does not show any clear variation with redshift. 
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1 INTRODUCTION 

Since the discovery of "ho les" in the X-ray emis sion at the cen- 
tre of the Perseus Cluster (Bohringer et al. 1993), and especially 
since the launch of Chandra and XMM-NEWTON, many depres- 
sions in the central intra-cluster medium (ICM ) of low redshift 
clusters have been found (e .g. Hydr a A, McNamara et al l200d : 



A2052 iBlanton et alj l200ll; A2 199. Ijohnstone et al.l 120021 : Cen- 
ta urus. [Sanders & Fabianl 1 20021). Recent compi latio ns are given 
bvlDunn & Fabianl feOOfj )- iRaffertv etafl feOOrj) and iBirzan etai] 
(2004). These holes have been observed to anti-correlate spectac- 
ularly with the radio emission from the active galactic nucleus 
(AGN) at the centres of these clusters. Their morphology, partic- 
ularly in the closest clusters, has led to their interpretation as bub- 
bles of relativistic gas blown by the AGN into the thermal ICM. 
The relativistic gas is less dense than the ICM, so the bubbles de- 
tach fro m the core and rise up buoyantly through the cluster, e.g. 
Perseus dChurazov etal]|20od IFabian et alj|2003l) . The older bub- 
bles tend not to have GHz radio emission associated with them and 
have been termed "Ghost" bubbles. 

The X-ray emission of the ICM represents an energy loss, 
so leads the plasma cooling in the core if there were no com- 
pensating energy source. To maintain pressure support, the gas 
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would flow on to the central galaxy as a "cooling flow." How- 
ever, with the high spatial and spectral resolution of Chandra and 
XMM-Newton, the gas temperature is seen to drop by a factor 
of several in the cor e bu t no continuous cooling flow is found 
dPeterson et alj [20031 see IPeterson & Fabian 2006 for a review). 
Some form of g entle, continuous an d distributed heating appears 
to be at work JVoigt & FabiarJl2004l) . The interaction of the AGN 
with the ICM is the likeliest explanation as the heat source for 
cool-core clusters. The bubble mo del for radio sources was pro- 
posed bv lGull & NorthoveJ fl973h . and further theoretical results 
were presented in Churazov et alj ( 1200 . 2002). Recent studies o n 
the heating ab i lity of AG N within clusters are feirzan et alj d2004l) : 
IRaffertv et ail bOOfj) and lDunn & Fabianl ( feOOfj) . 

A major ity (71 per c ent) of "cooling core" clusters harbour 
radio sources (Burns 1990), and a similar number of clusters which 
require heating (l i kely to be a cooling core) harbour clear bubbles 
dDunn et al. 2005: Du nn & FabianlioOrj) . The action of creating the 
bubbles at the centre of the cluster by the AGN is a favoured method 
of injecting energy into the central regions of the cluster and so 
prevent the ICM from cooling. This process sets up sound/pressure 
waves in the IC M, the dissipatio n of which requires the ICM is 
to be viscous ( IFabian et al1l2003l) . The viscous dissipation of the 
pressure waves allows the energy from the bubble creation to be 
dissipated far from the cluster centre where it is required. 

In lDunn & Fabianl 1200(3) we analysed clusters in the Bright- 
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est 55 (B55) sample tedge et alj[l990h . These are all at compara- 
tively low redshif t (only two are at z > 0.1). Using the results of 
IPeres et 

ID GUI) 

and the literature at the time we selected those 
clusters with both a short central cooling time and a large central 
temperature drop, as well as those with a central radio source. We 
found that at least 70 per cent of clusters in which high cooling 
rates are expected (those with short central cooling times and cen- 
tral temperature drops) harbour AGN bubbles. We also found that 
95 per cent of these cool core clusters harbour radio cores. The to- 
tal fraction of clusters which harbour a radio source is 53 per cent 
(29/55). Where AGN bubbles are present, on average they show 
that the AGN is injecting enough energy into the central regions 
of the cluster to offset the X-ray cooling out to the cooling radius 
(r C ooi where t coo i = 3 Gyr). To increase the number of clusters 
studied and also investigate any evolution in the AG N heating with 
redsh ift we now turn to the Brightest Cluster Sample ( lEbeling et al.l 
[l998j), initially performing a similar analysis and then combining 
the results from both cluster samples. 

The sample selection for this work is outlined in Section [2] 
and the data preparation and reduction in Section [3] We then look 
at the two different cluster subsets, those with bubbles in Section 
[4] and those with central radio sources in Section [5] The results of 
our investigations into the global differences within the different 
subsets is presented in Section|6] The combined cluster sample and 
the redshift evolution results are presented in Section [7] We use 
H = 70kms _1 Mpc -1 with fl m = 0.3, = 0.7 throughout 
this work. 



2 SAMPLE SELECTION 

To increase the number of clusters studied and to investigate the 
redshift evolution of AGN in cl usters we now turn to the ROSAT 
Brightest Cluster Sample (BCS, lEbeling et alj[l998l) . Most of the 
B55 clusters are at low redshifts (z < 0.1), whereas the BCS 
clusters extend to higher redshifts. The BCS is a 90 per cent flux- 
complete sample comprising the 201 X-ray-brightest clusters in 
the northern hemi sphere at high galactic latitudes (| b > 20° 
lEbeling etal.lll998h . The member clusters have measured redshifts 
of z < 0.4 and fluxes greater than 4.4 x 10~ 12 erg cm" 2 s" 1 in 
the 0.1 - 2.4 keV band (for Ho = 50 km s" 1 Mpc" 1 ). There is a 
low flux extension to the BCS, the eBCS, which includes clusters 
with fluxes in the range 2.8 x 10~* 2 to 4 .4 x 10~ 12 erg cm -2 s -1 
jEbeling et al j2000l) . Bauer et alj Jiool) investigated the existence 
of cool cores in the high redshift end of the BCS (z > 0.15), and 
found that 34 per cent of clusters exhibited signs of strong cooling 
(central t coo i < 2 Gyr). In our study (z > 0.1 and i coo i < 3 Gyr) 
we find that 16 clusters have a short central cooling time, corre- 
sponding to 38 per cent. If the temperature drop is also required 
for "strong cooling" then the fraction drops to 33 per cent, both of 
which match their result. Globally there have been a v ariety of esti- 
mates on the fraction of clusters which have cool corefl Per es et al.l 
d 19981) used the B55 sample and found a fraction of 70-90 per 
cent. However they define a "cooling flow" cluster is that which 
has t coo i < ^duster = 13 Gyr. If we use t coo \ < 13 Gyr for our 
selection of cool core clusters, 38/42 (90 per cent) of clusters have a 
cool core. It is therefore reasonable to assume that the Chandra ob- 
servations of clusters within the BCS have not been highly biased 

1 The ones which would previously have been described as having a "cool- 
ing flow". 



toward s cool core clusters. In the B55 sample from lDunn & Fabianl 
(2006) only 1/30 clusters would not have a cool core, though the 
selection procedure was different in that study. 

So that the BCS and B55 cluster samples do not overlap in this 
study, we only select those clusters in the BCS with z > 0.1 as the 
parent sample for this study. There is only one cluster in common 
with both the B55 and BCS with this cutoff, A2204. There are 87 
clusters (including A2204) in the parent sample, of which 42 have 
observations in the Chandra archive. These are shown in Table [T] 
We assume that the Chandra observations have not been biased to 
a particular type of cluste^f] and so quote population fractions from 
this sample of 42 rather than from the parent sample of 87 clusters. 

Following the data reduction (see Section[3]l this parent sam- 
ple wa s further split into sub-samples, following iDunn & Fabianl 
(2006): those clusters which harbour clear bubbles; those in which 
heating is required to prevent rapid cooling and those which har- 
bour a radio source were identified. The allocation of clusters into 
the three sub-samples is shown in Table[T] 

To identify those clusters which are likely to require some 
form of heating to prevent large quantities of gas from dropping out 
we take those which have a short central cooling time (< 3 Gyr) 
and a large central temperature drop (T c<m tre/7outer < 1/2). Ex- 
actly half (21) of the clusters in the sample have a central temper- 
ature drop, an d 16 have a short c entral cooling time (16/42, 38 per 
cent, similar to lBauer et al.l2005l although they use t coo i < 2 Gyr). 
14 clusters have both (14/42, 33 per cent), of which at least 6 
have clear bubbles (6/14, 43 per cent; 6/42, 14 per cent). We use 
the NVSSEI to determine whether the clusters have a central ra- 
dio source. 23 clusters harbour a central radio source (23/42, 55 
per cent), including all those which require some form of heating. 
There appears to be a category of clusters which have cool cores 
and a radio source but in which no clea r bubbles have been de - 
tected. Both in the B55 sample (Table 2 IDunn & Fabianl d2006l) . 
2 column) and the BCS have these sorts of clusters, with similar 
fractions in each (5/30 in the B55 and 8/4 2 in the BCS). 

Comparing with the B55 sample, IDunn & Fa bian fe006l) 
found that the fraction of clusters which require heating is 25 per 
cent, and which harbour a radio source is at least 53 per cent. Of 
the clusters which require some form of heating, at least 70 per- 
cent harbour clear bubbles and 95 percent harbour a central radio 
source. 

In comparing the number of identified bubbles in the BCS and 
B55 sample it is expected that as a result of the greater distance 
of clusters in the BCS, fewer bubbles would be identified even if 
the same proportion of clusters harboured them. It is therefore not 
surprising that the percentage of clusters which harbour bubbles is 
half that determined for the B55 (from those clusters which require 
heating). The proportion of clusters which harbour radio sources, 
on the other hand, is very similar (BCS=55, B55=53 per cent) as is 
the proportion of clusters requiring heating (BCS=33, B55=25 per 
cent). 



2 Only around 1/3 of the clusters observed with Chandra would be classed 
as cool-core clusters. A quick analysis of the average length of observations 
for different cluster types is discussed in Section[6] 

3 NRAO (National Radio Astronomy Observatory) VLA (Very Large Ar- 
ray) Sky Survey 
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Table 1. Cluster Sample 



Cluster 


Redshift 


Obsld Exposure (ks) 




T drop 


Both 


Bubbles 


Radio 


A68 


0. 


255 


3250 


9.9 












A115N 


0. 


197 


3233 


49.7 


Y 


Y 


Y 


Y 


Y 


A115 S 


0. 


197 


3233 


49.7 




Y 








A267 


0. 


227 


1448 


7.4 












A520 


0. 


202 


528 


9.5 




Y 








A586 


0. 


171 


530 


10.0 










Y 


A665 


0. 


182 


531 


8.9 












A697 


0. 


282 


4217 


19.5 












A773 


0. 


217 


5006 


19.8 












A781 


0. 


298 


534 


9.9 




Y 






Y 


A963 


0. 


206 


903 


35.9 




Y 








A1068 


0. 


139 


1652 


25.5 


Y 


Y 


Y 




Y 


A1201 


0. 


169 


4216 


24.0 












A1204 


0. 


171 


2205 


23.3 


Y 










A1413 


0. 


143 


5003 


75.0 












A1423 


0. 


213 


538 


9.7 


Y 


Y 


Y 




Y 


A1682 


0. 


226 


3244 


5.4 




Y 






Y 


A1758 N 


0. 


279 


2213 


49.7 










Y 


A1763 


0. 


223 


3591 


19.5 










Y 


A1835 


0. 


252 


495 


19.3 


Y 


Y 


Y 


Y 


Y 


A1914 


0. 


171 


3593 


18.8 












A2111 


0. 


229 


544 


10.2 




Y 








A2204 


0. 


152 


499 


10.1 


Y 


Y 


Y 




Y 


A2218 


0. 


171 


1666 


40.1 












A2219 


0. 


228 


896 


42.3 










Y 


A2259 


0. 


164 


3245 


10.0 












A2261 


0. 


224 


5007 


24.3 


Y 










A2294 


0. 


178 


3246 


9.0 




Y 








A2390 


0. 


231 


4193 


89.8 


Y 


Y 


Y 




Y 


A2409 


0. 


147 


3247 


10.2 












HercA 


0. 


155 


6257 


49.4 


Y 


Y 


Y 


Y 


Y 


MS0906+1110 


0. 


175 


924 


29.6 












MS 1455+2232 


0. 


258 


4192 


91.4 


Y 


Y 


Y 




Y 


RXJ0439+0520 


0. 


208 


1649 


9.6 


Y 


Y 


Y 




Y 


RXJ0439+0715 


0. 


245 


3583 


19.2 










Y 


RXJ1532+3021 


0. 


362 


1649 


9.2 


Y 


Y 


Y 


Y 


Y 


RXJ1720+2638 


0. 


161 


4361 


23.3 


Y 


Y 


Y 




Y 


RXJ2129+0005 


0. 


234 


552 


9.9 


Y 


Y 


Y 




Y 


ZwC11953 


0. 


373 


1659 


20.7 










Y 


ZwC12701 


0. 


214 


3195 


26.7 


Y 


Y 


Y 


Y 


Y 


ZwC13146 


0. 


285 


909 


45.6 


Y 


Y 


Y 


Y 


Y 


ZwC15247 


0. 


229 


539 


9.0 










Y 


Totals 








42 


16 


21 


14 


6 


23 



All the clusters in the sample, for the sample selection see text. The exposure time is that after reprocessing the data. The cooling times and central temperature 
drops have been determined from the profiles calculated during the course of this work. Clusters with < 3 Gyr are classed as having a short central cooling 
time, and those with T ce ntre/7outer < 1/2, a central temperature drop. The radio detections come from the NVSS. 



3 DATA PREPARATION 

The X-ray Chandra data of the clusters were processed and cleaned 
using the CIAO software and calibration files (CIAO v3.3, CALDB 
v3.2). We began the reprocessing by removing the afterglow detec- 
tion and re-identifying the hot pixels and cosmic ray afterglows, 
followed by the tool ACIS_PROCESS.EVENTS to remove the pixel 
randomisation and to flag potential background events for data ob- 
served in Very Faint (VF) mode. The Charge-Transfer Efficiency 
was corrected for, followed by standard grade selection. Point- 
sources were identified using the WAVDETECT wavelet-transform 
procedure. For clusters observed with the ACIS-S3 chip, the ACIS- 
Sl chip was used to form the light curves where possible. In all 
other cases, light-curves were taken from on-chip regions as free 



as possible from cluster emission. For the spectral analysis, back- 
grounds were taken from the CALDB blank-field data-sets. They 
had the same reprocessing applied, and were reprojected to the cor- 
rect orientation. 

Cluster centroids were chosen to lie on the X-ray surface 
brightness peak. Annular regions, centred on the centroids, were 
automatically assigned with constant signal-to-noise, stopping 
where the background-subtracted surface brightness of the cluster 
dropped below zero. The initial signal-to-noise was 100, and this 
was increased or decreased by successive factors of v2 to obtain a 
number of regions between four and ten. The minimum signal-to- 
noise allowed was 10. For clusters imaged on the ACIS-I chip, the 
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Figure 1. The distribution of i coo i in the BCS sample. The black bars in- 
dicate those clusters which harbour clear bubbles, and the grey ones, those 
with a central radio source. ZwC15247, which harbours a radio source, is 
beyond the right edge of this plot. 



effects of excluding the chip-gaps was accounted for by adjusting 
the stop angles of the annuli. 

The 0.5 — 7keV spectra were extracted, binned with a mini 



mum of 20 counts per bin, and, using XSPEC (vl2.3) ("e.g.lArnauc 
1996), a PROJCT single temperature MEKAL (e.g. iMewe et al 



1995) model with a PHABS absorption was used to deproject the 
cluster. Chi-square was used as an estimator of the spectral fit. In 
some clusters the temperatures for some of the regions were unde- 
fined. To solve this problem, the minimum number of regions was 
reduced, a maximum radius for the outermost annulus was set, or 
the outermost region was removed to try to improve the behaviour 
of the profile. Using the deprojected cluster temperature, abundance 
and normalisation profiles; density, pressure, entropy, cooling time 
and heating profiles were calculated. 

These profiles give azimuthally averaged values for the cluster 
properties and have been used in the subsequent calculations. In 
some clusters, for example M87 and Perseus, the central parts of 
the clu ster are not very smooth, e.g. due to bubbles. [Donahue et al.l 
( 2006) show that these features do not strongly bias estimates of the 
entropy. As such the use of these azimuthally averaged values is not 
likely to introduce large biases into the subsequent calculations. 



4 CLUSTERS WITH BUBBLES 

The analysis of the cl usters with clear bubbles is the same as in 
iDunn & Fabianl J2006l) . Where clear cavities are seen in the X-ray 
emission then the energy injected into the central regions of the 
cluster is calculated from the energy in the bubbles. 



-Ebut 



7i 



PthV + pthV 



7i 



7i 



-pthV 



(1) 



where V is the volume of the bubble, p t h is the thermal pressure of 
the ICM at the radius of the bubble from the cluster centre and 71 
is the mean adiabatic index of the fluid in the bubble. Assuming the 
fluid is relativistic, then 71 =4/3, and so -Ebubbie = 4p t h V. 

We parametrise the bubbles as ellipsoids, with a semi-major 
axis along the jet direction, Ri, and a semi-major axis across it, 
R w . The depth of the bubbles along the line of sight is unknown 



and so we assume them to be prolate ellipsoids, with the resulting 
volume being V = 4ttR^,Ri/3. 

To estimate the rate at which the AGN is injecting energy into 
the cluster the we need to estimate the age of the bubble. There 
are two principle methods of estimating the age of the bubble, the 
sound speed timescale and the buoyancy rise time. Observations of 
nearby clusters have shown that young bubbles which are still being 
inflated do not create strong shocks during their inflation. They are 
therefore assumed to be expanding at around the sound speed, and 
their age is t agc — tc s = Ri/c s where the sound speed in the ICM 
is 



c s = \/72£;b T/fimu 



(2) 



where 72 = 5/3. As bubbles evolve they are expected to de- 
tach from the AGN and rise up buoyantly through the ICM 
dChurazov et alj f2000). Their age is best calculated by the buoy- 
ancy rise time; t agc — £buoy = iidist / «b UO y , with the buoyancy 
velocity given by 



v huoy = y/2gV/SCr> (3) 

where Cd = 0.75 is the drag coefficient dChurazov etai]|200lh 
and S = ttR^, is the cross-sectional area of the bubble. 

The radio sources in this cluster sample are analysed after the 
identification of cavities within the X-ray emission. We find two 
strong radio sources within this sample, A115 (3C 28) and Her- 
cules A (3C 348). Both have current radio emission in a bi-lobed 
morphology, and the radio lobes are still close to the core. None of 
the others have radio emission correlated with the cavities in the X- 
ray emission. We therefore assume that the radio emission has aged 
and so faded beyond the sensitivity limit, and class these bubbles 
as Ghost bubbles. 

Using the deprojected temperature, density and abundance 
profiles, X-ray cooling was calculated for each spherical shell. This 
was compared to the power injected into the cluster by the AGN 
in the form of bubbles (Pbubbic = Eb u bbie/ta g e)- Tne distance 
out to which the bubble power can counteract the X-ray cooling 
was therefore calculated. Uncertainties have been estimated using 
a simple Monte Carlo simulation of the calculatioiQ 

As this cluster sample is at a higher average red shift than the 
Brightest 55 sample studied in lDunn & Fabianl |2006), fewer clear 
bubbles have been identified. Also, the X-ray emission is fainter, 
and so analysis of the outermost regions of the cluster becomes 
more difficult. Only large or extreme bubbles will be detected. It is 
therefore not surprising that of the six clusters which have clear de- 
pressions in the X-ray emission, in three of them the AGN is inject- 
ing much more energ y than is being los t through X-ray emission in 
the central regions. iNulsen et alj d2005l) have studied the Chandra 
data of Hercules A and concluded that this outburst is exception- 
ally powerful. A comparison between the X-ray cooling luminosity 
within the cooling radius to the bubble power is shown in Table[3] 

The average distance, as a fraction of the cooling radius, out 
to which the bubbles can offset the X-ray cooling calculated from 
the results is presented in Table [3] The uncertainty in the mean 
was estimated using a simple bootstrapping method. On average, 
the clear bubbles in this sample counteract the X-ray cooling out to 

4 For this and the calculation of the expected bubble sizes in Section [3] 
the method is as follows. The input values were assumed to have a Gaus- 
sian distribution. The values used in each run of the simulation were then 
randomly selected. After a large number of runs (2 X 10 5 ) the resulting pa- 
rameters were sorted into ascending order and the median and interquartile 
range are used as the value and uncertainties. 
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Table 3. Clusters with Bubbles: AGN Heating ability 



Cluster 


L a 

cool 


4PthV/tagc 


^"cool 


f*hcat 


r hcat / r CO ol 


Heated Fraction 2 




(10 43 ergs" 1 ) 


(10 43 ergs" 1 ) 


(kpc) 


(kpc) 






A115N 


16.8 ± 1.1 


273 ± 33 


70.5 ± 3.8 


> 153 


2.17 ± 0.12 


16.4 ± 2.2 


A1835 


205 ± 6 


190 ± 38 


96.7 ± 2.3 


90.5 ± 14.9 


0.94 ± 0.16 


0.92 ± 0.19 


HercA 


4.2 ± 0.4 


1560 ± 360 


45.9 ± 3.4 


> 600 


13.1 ± 1.0 


373 ± 98 


RXJ1532+30 


218 ± 23 


331 ± 126 


107 ± 10 


187 ± 96 


1.75 ± 0.91 


1.52 ± 0.60 


ZwC12701 


23.4 ± 1.9 


602 ± 82 


67.4 ± 3.0 


> 301 


4.47 ± 0.20 


25.5 ± 4.1 


ZwC13146 


153 ± 7 


786 ± 144 


93.0 ± 3.7 


> 500 


5.38 ± 0.21 


5.11 ± 0.97 



The clusters in bold only have a lower limit for r^eat /r coo \ as the AGN counteracts all the cooling from the X-ray emission which could be analysed. a 
L coo l is calculated for the 0.5 — 7.0 keV range. b The radius rjjeat. as a fraction of the cooling radius r coo i, out to which the energy injection rate of the 
bubbles (4p t hV /4 a gc) can offset the X-ray cooling within that radius. c The fraction of the X-ray cooling that occurs within r[ lcat . 



Table 2. Clusters with Bubbles: Bubble Properties 



Cluster 


Bubble 


Type 


Ri 


R w 


^dist 








(kpc) 


(kpc) 


(kpc) 


A115 


NE 


Y 


32.1 ±3.2 


28.78 ±2.8 


39.0 


±3.9 




SW 




27.2 ±2.7 


20.3 ±2.0 


45.2 


±4.5 


A1835 


NE 


G 


15.5 ± 1.6 


11.6 ±1.2 


23.3 


±2.3 




SW 




13.6 ±1.4 


9.6 ± 1.0 


16.6 


±1.7 


HercA 


E 


Y 


110 ±11 


75.4 ±7.5 


110 


±11 




W 




124 ± 12 


75.4 ±7.5 


124 


±12 


RXJ1532+30 


W 


G 


23.5 ±2.4 


23.5 ± 2.4 


32.0 


±3.2 


ZwC12701 


E 


G 


29.8 ±3.0 


29.8 ±3.0 


42.6 


±4.3 




W 




28.7 ±2.9 


28.7 ±2.9 


39.5 


±4.0 


ZwC13146 


NE 


G 


28.0 ±2.8 


28.0 ±2.8 


57.1 


±5.7 




SW 




30.0 ±3.0 


30.0 ±3.0 


53.2 


±5.3 



The bubbles are parametrised as ellipsoids, with R\ the semi-major axis 
along the jet direction, and R v the semi-major axis across it. We class 
young bubbles as those with clear radio emission, and ghost bubbles as 
those without. Both A115 (3C 28) and Hercules A (3C 348) have strong 
bi-lobed radio sources at their centres. 



''heat /''cool = 1.34±0.40. The fraction of the X-ray cooling within 
the cooling radius offset by the action of the AGN is 1.22 ± 0.15. 
We emphasise that this average is only for two clusters (A1835 and 
RXJ1532+30), and should be regarded as a lower limit for the fol- 
lowing reasons. 

The values exclude those clusters where the calculation only 
provides a lower limit on the distance. Adding these clusters into 
the calculation results in ri lca t /r coo \ = 4.31 ± 1.10 and a heated 
fraction of 75.0 ± 50.4. It is likely that, because of the greater av- 
erage distance of this cluster sample, that the clusters with clear 
bubbles are those which are undergoing/have recently undergone 
powerful AGN outbursts. The more gentle energy injection events 
would not be easily identified in the X-ray emission. 

We do not comment further on these results here, as there are 
only six clusters in this sub-sample. Further analysis, where these 
clusters are combined with those from the B55 sample, is discussed 
in Section|7] 



5 CLUSTERS WITH RADIO SOURCES 

All the clusters which require some form of heatin g harbour a radio 



source , which is similar to the fraction found by 
20061) of 95 per cent. Following the study in 



Dunn & Fabian 



Dunn & Fabiar 



2006) we invert the problem, and use the results from the clus- 
ters with bubbles on this sample. Using the average distance out to 
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Figure 2. The 8 GHz image of A 1763 from the VLA in A-array. This is 
the only cluster in which bubbles are expected but none seen in the X-ray 
emission, which also has a bi-lobed radio source at the centre where the 
structure of the radio emission was resolved in the archival VLA data. 



which the bubbles counteract the X-ray cooling, we calculate the 
X-ray cooling luminosity for these clusters and equate this to the 
energy injected by the AGN. 

-Ftubblc = L coo \(r < Theat), 

where rh ca t is calculated from rhcat = 1.34r CO oi. We assume that 
these "expected" bubbles are young, and so are being inflated at the 
sound speed. 



Rb 



4pthC s 



Uncertainties have been estimated using a simple Monte Carlo sim- 
ulation of the calculation. 

Using the size of an expected bubble in these clusters, the X- 
ray count total within the bubble was calculated from the current 
X-ray observation. This allowed the signal-to-noise of the observa- 
tion to be obtained The difference in counts between the bubble 

5 The signal-to-noise is calculated from the predicted size of the bubbles 
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Table 4. Predicted Bubble Sizes 



Cluster 


z 


f*cool 


Predicted Bubble Radius 


Central 


Radio dimensions 






(kpc) 


(kpc) 


(arcsec) 


X-ray S/N 


(arcsec) 


A586 


0.171 


20.59 ± 10.63 


2.39 ± 1.00 


0.82 ±0.34 


1.7 




A781 


0.298 


2D Qfi + 8 23 


9.75 ± 2.04 


2.18 ±0.46 


1.3 


5.2 X 5.2 


A1068" 


0.139 


77.19 ± 1.86 


14.51 ± 2.23 


5.91 ±0.91 


59.6 


5.5 x 5.5 


A 1423 


0.213 


46 02 + 8 68 


5.89 ± 1.69 


1.69 ±0.49 


3.1 


3 9x103 


A1682 


0.226 


27.82 ± 8.79 


5.87 ± 1.07 


1.61 ±0.29 


1.1 


5.8 x 5.8 


A1758 N 


0.279 


12 37 + 12 17 


4.82 ± 1.68 


1.13 ±0.40 


2.6 




A 1763 


0.223 


1 8 38 + 4 24 


4 ™ + n w 

rt.vJi? _l_ \J.tJiJ 


i 22 + 1 5 


2.0 


1.4 x 3.4 


A2204 


0.152 


69.47 ±3.19 


9.08 ± 3.57 


3.42 ±1.34 


38.0 


5.0 x 5.0 


A2219 


0.228 


20.55 ± 2.34 


4.64 ± 0.29 


1.26 ±0.08 


5.3 


4.6 x 4.6 


A2390 


0.231 


60.91 ± 1.16 


10.81 ± 2.64 


2.92 ±0.71 


30.8 




MS 1455+2232 


0.258 


96.56 ± 1.10 


19.62 ±4.26 


4.88 ±1.06 


65.4 




RXJ0439+0520 


0.208 


74.66 ± 6.20 


11.81 ±3.41 


3.45 ±1.00 


13.8 




RXJ0439+0715 


0.245 


40.18 ± 3.51 


5.45 ± 0.29 


1.41 ±0.08 


3.3 




RXJ1720+2638 


0.161 


79.42 ± 2.77 


12.64 ± 1.36 


4.54 ±0.49 


32.1 


amorphous source 


RXJ2 129+0005 


0.234 


72.78 ± 5.00 


9.25 ± 2.86 


2.47 ±0.76 


8.5 




ZwC11953 


0.373 


22.06 ± 8.33 


4.48 ± 0.94 


0.86 ±0.18 


1.7 


3.7 x 3.7 


ZwC15247 6 


0.229 


< 78.54 


< 8.83 


< 2.4 


0.6 





a The radio emission from A1068 is discussed in more detail in the text. b The observation of ZwC15247 is short for such a distant cluster. It also does not 
have a cooling core, although it appears to have a central radio source, and so the X-ray emission is not strongly peaked at the centre. The deprojection analysis 
results in a cooling time which is not well constrained at the centre. We therefore give only the upper limits for the values obtained by the analysis for the 
quantities in this table. 



centre and the rims surrounding it was estimated from the bubbles 
in A2052 and Hydra A. The counts in the bubbles from these two 
clusters are around 30 per cent lower than the counts in the rims. If 
bubbles are present, then they would need to be detected above the 
noise at more than the 3<r level. For this to be the case, the signal- 
to-noise has to be greater than 13 (equivalent to noise at a 8 per cent 
level). However, as the cluster emission is changing over the scale 
of the bubbles the we may need a higher signal-to-noise to be able 
to detect any depressions. For a 5<r detection, the signal-to-noise 
required is 36. This does require knowing where the bubbles are so 
that the counts within the region where a bubble is expected can be 
compared to that outside. The results of this analysis are shown in 
Table E] 

There are six clusters (A1068, A2204, A2390, MS 1455+2232, 
RXJ01439+0520 and RXJ1720+2638) where the central signal-to- 
noise of the observation is such that if there are bubbles present 
within this region, then with the aid of any radio emission to guide 
the eye, the bubbles should be detectable at the 3a level. There 
are only three if the detection threshhold is lifted to 5a (A1068, 
A2204 & M S 1455), of which A106 8 is discussed below, A2204 is 
discussed in iDunn & Fabianl feOOfih and no high resolution radio 
data was found for MS 1455. As this signal-to-noise calculation has 
been calculated for the entire central part of the cluster, the value 
obtained may over estimate the signal-to-noise present at the loca- 
tion of the bubble if the core of the cluster is very bright. 

The bi-lobed radio emission A1068 appears offset to the X- 
ray cluster centroid. It would be expected that for a typical double 
source, the lobes would sit either side of the X-ray peak (under 
the assumption that the AGN is in the very centre of the cluster). 
However, one "lobe" coincides with the peak in the ICM emis- 
sion, the other occuring at 38 kpc from the centre. In the vicin- 
ity of this more distant lobe, the signal to noise is ~ 41. As 



and the background subtracted counts per pixel within the central annular 
region. 



no clear features are seen in the X-ray emission, and the lobe 
locations are not typical, we investigated what this radio source 
is coincident with. Using the Sloan Digital Sky Survey (SDSS, 
lAdelman-McCarthv & for the SPSS Collaboratio j 120071) we ex- 
tracted images of the sky centred on the brightest cluster galaxy 
(BCG). Aligning the radio and optical images showed that the more 
distant (south-western) lobe is coincident with another galaxy in the 
cluster (z = 0.136 compared to z = 0.138 for the BCG). The cen- 
tral "lobe" is therefore an unresolved radio core. This would also 
explain the lack of bubble features in the high signal-to-noise X-ray 
emission, as no bubbles are seen, given the lack of extended radio 
emission. 

Using the Very Large Array archive we attempted to 

find observations of the clusters in which bubbles are expected to 
investigate the morphology of the central radio source. We analysed 
data at a variety of frequencies and array configurations, depending 
what was available in the archive and which had the longest time- 
on-source. 

These clusters are all at redshifts higher than those in the B55 
sample, and so we expected to resolve fewer of the central radio 
sources. Table|4]shows the expected bubble sizes and measured ra- 
dio source dimensions. If no radio source dimensions are given then 
the source was not resolved in the observations that were analysed. 

The predicted bubble and observed radio source sizes are very 
similar in A1068, A1763 (Fig. Hi and A2204. A 1 068 is discussed 
earlier in this section and A2204 was discussed in lDunn & Fabianl 
J2006h and will not be covered further here. However, the radio 
source size and the X-ray observation signal-to-noise in A 1763 are 
such that detecting any X-ray features of the bubble would be im- 
possible, even with a much longer exposure. Under the assumption 
that the natural course of events when an AGN ejects radio emitting 



e The National Radio Astronomy Observatory is operated by Associated 
Universities, Inc., under cooperative agreement with the National Science 
Foundation. 
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jets into the central regions of a cl uster is that bubb l es are formed, 
then, as in the B55 cluster sample f5unn & Fabianl J2006I) . we are 
missing a large number of bubbles. This may be a timescale effect, 
currently in some of these clusters the AGN is not inflating bub- 
bles, or they are present but too small to be detected with either the 
current X-ray or radio data. 

In nearby clusters where X-ray and radio spatial resolution do 
not pose a problem for studying the morphology of the AGN and 
cluster, whenever there is a lobed radio source, there are always 
signs of an interaction between the relativistic and thermal plas- 
mas. Therefore, for more distant objects, where even the spatial 
resolution of Chandra is insufficient to allow the detailed investi- 
gation of the interaction, perhaps the radio emission alone will be 
enough to determine the action of the AGN on the cluster. If the 
two plasmas do not significantly mix over a timescale of ~ 10 s yr, 
then if the AGN is actively producing jets, bubbles would be the 
natural consequence, along with the associated heating of the ICM. 
Therefore, if the bubbles do allow a good calibration of the AGN 
energy injection rate, the dimensions of the radio emission rather 
than of any X-ray features could be used to calculate the heating 
from the AGN. 



6 CLUSTERS REQUIRING HEATING 

All the clusters which require heating (but do not have clear bub- 
bles - either reported or determined from the Chandra data) har- 
bour a central radio source (as determined from the NVSS). Their 
expected bubble sizes have been calculated in Section[5] 

As a comparison between the three different sub-samples, Fig. 
[3] shows the temperature, entropy and cooling time profiles for the 
different sub-samples. There are only six clusters in the sub-sample 
containing clear bubbles and so small number statistics are relevant. 
Fig. [3] also shows the best fitting powerlaw slope to all the data 
points. It is interesting to note that in the average slope for all three 
parameters, the slope for the clusters with radio emission is always 
in between those with bubbles and those with neither. The temper- 
ature slopes of the radio and no-radio clusters are almost identical, 
whereas with the entropy and cooling time, it is the bubbled and 
radio clusters which have very similar slopes. 

The clusters which require some form of heating are high- 
lighted in Fig. [3] The clusters which require some form of heat- 
ing do appear to be the coolest clusters and those which have the 
lowest central entropy. They are naturally those with the shortest 
central cooling time as this was one of the selection criteria. As 
is clear from Fig. [3] these clusters are those in which t he deprojec- 
tion re sulted in the smallest central bin. As discussed in lBauer et al.l 
(2005), the radius of the innermost annulus depends on the cluster 
distance and the signal-to-noise of the observation. The latter pa- 
rameter is related to whether the cluster exhibits a cool core or not, 
as cool core clusters have higher central densities and hence have 
brighter X-ray cores. 

There is a further possible bias towards having detailed/deep 
observations of cool core clusters and less deep observations of 
other clusters. A quick calculation of the average exposure for clus- 
ters which have bubbles, a radio source or neither results in 33.5, 
28.4 and 24.3 ks respectively. Although these exposure times are 
not significantly different (standard deviation ^ 17) it is interest- 
ing to note that the clusters which are of most interest to the heat- 
ing and cooling balance are those which appear to have the longest 
exposure times. These clusters, as they are cool core, are also the 
brightest. For a balanced analysis the non-cool core clusters prob- 



ably should have a longer average exposure times as they are less 
bright. 

If the entropy data are split into those clusters which require 
some form of heating and those that do not, then there is a large 
difference in the slopes of the average entropy profiles. The dif- 
ferent best fitting slopes are 1.01 and 0.72 respectively. Ignoring 
A 1204 in the entropy fits for the clusters with no radio data results 
in a slope of 0.77. A1204 only just fails on the temperature drop 
selection criterion for requiring heating. 

The average cooling time profiles of the clusters with bub- 
bles and those with a radio source match the profiles presented in 
IVoigt & Fabianl fefXMl) of tmn\ oc r 1 3 . 



7 DISCUSSION 

7.1 Combining the Samples 

In order to study a large number of clusters over a range of redshifts 
we combine the BCS and the B55 samples. As A2204 is present in 
both we have used the central cooling time from the analysis pre- 
sented in this manuscript in the following discussion. We therefore 
have 30 + 42 — 1 = 71 clusters in the combined sample. 

The B55 and BCS samples have slightly different selection 
criteria. The BCS sample selects clusters with S > 0°, \b\ ^ 
20°, z < 0.3 and S0.1-2.4kcV > 4.4 x 10~ 12 erg cm -2 s _1 , 
with the extended sample down to So.i-2.4keV > 2.8 x 
10 _12 ergcm -2 s _1 . The B55 sample has S2-iok c v > 1.7 x 
10~ n erg cm -2 s _1 as its criteria and includes 9 clusters at low 
galactic latitudes (|6| < 20°). Both are around 90 per cent flux 
complete. We do not adjust the selection criteria for consistency 
between the two samples. 

So that both the BCS and B55 cluster samples have the same 
data reduction and spectral fitting applied we have re-reduced the 
B55 cluster data. This has allowed us to take advantage of the im- 
proved data-reduction script - e.g. taking into account the chip gaps 
in ACIS-I observations. 



7.7.7 Cooling times 

Fig.|4]shows the distribution of cluster central cooling times. Rather 
than t ake the values from the analysis of ROSAT data lPeres et al.l 
( 1998), we have repeated the extraction of the central cooling time 
performed here on the BCS sample, on the B55 sample. We note 
that not all of the clusters in the high redsh ift end of the BCS sam- 
ple have been observed with Chandra; in iDunn & Fabianl (2006) 
clusters without short central cooling times or central radio sources 
were discarded from further analysis, whereas all clusters for which 
Chandra observations exist have been analysed in this work. As a 
result the fraction of clusters harbouring bubbles/radio sources as a 
function of the central cooling time are estimates. 

What can be seen from Fig. [4] however, is that radio sources 
(with or without clear bubbles) are much more common in clusters 
with central cooling times of 4 Gyr or less. In fact there are only 
10/24 clusters with i coo i ^ 4 Gyr which have central radio sources. 
In the distribution with finer binning, the peak in the central cooling 
times appears at around 0.6 Gyr. In the current combined samples, 
there appear to be no clusters which harbour clear bubbles and have 
central cooling times > 1.2 Gyr. Very few clusters (5/47) with cen- 
tral cooling times shorter than 4 Gyr have no evidence for a central 
radio source (from the NVSS). 
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Figure 3. LEFT: Temperature, CENTRE: Entropy and RIGHT: Cooling time profiles for TOP: Clusters with clear bubbles, MIDDLE: Clusters with central radio 
sources and BOTTOM: clusters with no central radio source. Best fitting powerlaw slopes, (r), to all the date on the plots are shown. The clusters which require 
heating are highlighted by a large circle in the plots for the clusters harbouring a central radio source. 



Although the two samples combined here have not been cre- 
ated in the same way (lack of non-cool core/non-radio source har- 
bouring clusters in the B55 analysis) there appears to be an excess 
of clusters which have a cooling time of < 1 Gyr. To take into ac- 
count the effect of all the clusters in the B55 sample which have not 
yet been analysed on Fig [4] we look at the distribution of t coo \ at 
> 3 Gyr. The distribution is fairly flat, so even if all the remaining 
clusters follow this fiat distribution, then the "pile-up" of clusters 
at tcooi < 1.2 Gyr would still be seen. 

An explanation of this pile-up could result from the feedback 
cycle suggested as a solution to the cooling flow problem. Only 
once the central cooling time is less than ~ 1 Gyr will the AGN 
be active and inject energy into the central regions of the cluster 
and create the tell-tale bubbles. It would depend on whether the 
AGN activity is gentle and relatively continuous or explosive and 



intermittent as to what the final state of the cluster would be. An ex- 
plosive and energetic AGN outburst would increase the ICM tem- 
perature, and hence cooling time by a large amount. As a result 
the next injection event would occur after a long dormant period. 
However, a gentle injection of energy would only reheat the ICM 
a small amount, and so the next injection event would take place 
relatively shortly thereafter. 

The large number of clusters at short cooling times, some of 
which have radio sources but no clear bubbles, implies that the ma- 
jority of injection events are the "gentle and often" kind; as most 
clusters with radio sources (dormant AGN?) and bubbles are found 
with short central cooling times. 

iBauer et al] ( l2005l) investigated the cooling time profiles of the 
Chandra-obseived BCS clusters. There are five more clusters in our 
sample as we use a lower minimum redshift for our sample selec- 
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Figure 4. TOP: The distribution of t coo i in the combined cluster sample. 
BOTTOM: A finer binning of the clusters with a short central t coo i. The 
black bars indicate those clusters which harbour clear bubbles, and the grey 
ones, those with a central radio source. 

tion and there have been further observations of BCS clusters. We 
compare our central cooling times to theirs and find that on aver- 
age ours are shorter by around a factor 0.7. There are some differ- 
ences between the data-reduction methods of the two studies. We 
d o use the cluster ce ntroid as calculated by the CIAO software as 
in lBaueretalJj2005» . but then visually inspect the annular regions 
and manually change their location to match the surface brightness 
peak (12 clusters). We also allow the galactic abso rption to be un- 
constrained in the spectral fitting. iBauer et al] d2005t) state that us- 
ing only the centroid rather than the surface brightness peak could 
bias their central cooling times up by up to a factor 3. 

7.7.2 Heating-Cooling balance 

We add the six clusters in the BCS sample which harbour clear 
bubbles to those 16 f rom the B55 sample and update Fig. 2 of 
Dunn & Fabian ( 2006) for the cooling luminosity and the bubble 
power (Fig. [5). The colour scale shows the redshift of the cluster. 
The six new clusters are all at the top of the plot in red as they are 
all at z > 0.15. 

The average distance, as a fraction of the cooling radius, out 
to which the bubbles can offset the X-ray cooling from the com- 




£ C ooi(erg/s) 



Figure 5. The rate of energy loss via X-ray emission from within the cool- 
ing radius of the cluster versus the bubble power from the AGN. These are 
for the clusters which harbour clear bubbles. The two triangle symbols are 
for M87 and Perseus where L coo i could not be determined from a single 
ACIS chip. The colour scale is the logarithm of the redshift, blue is low 
redshift, red is high redshift. The line is for exact balance between Pbubblc 
and L cool . 



bined sample is rteat /Vcool = 0.92 ± 0.11. The fraction of the 
X-ray cooling within the cooling radius offset by the action of the 
AGN is 0.90 ± 0.13. We ignore those clusters where the bubbles 
are sufficiently powerful to offset all the cooling within the anal- 
ysed region, and so only a lower limit on the distance is obtained. 
Including these clusters gives r^eat /Vcool = 2.10 ± 0.40 and the 
fraction of 20.92 ± 15.11. Although there are selection effects, es- 
pecially at the high redshift end of the sample, we can conclude 
that on average the bubbles counteract the X-ray cooling within the 
cooling radius. 

The problem with this plot is that it shows an instantaneous 
measure of L coo i and Pbubbie rather than a time average. However 
for the lower redshift sources it is still reasonable to say that on av- 
erage the bubble power is sufficient to counteract the X-ray cooling. 
At the higher redshifts all clusters lie on or above the line of exact 
balance. Under the assumption that clusters at redshifts of around 
0.2 — 0.3 are similar to those in the local Universe, we would ex- 
pect to see some clusters which have bubbles which appear to be 
insufficient to counteract all the X-ray cooling. 

Fig. [6] shows that there is a trend between the average size of 
a bubble and the X-ray cooling luminosity. The smaller L coo i, the 
smaller the detected bubbles. This trend does point to some form of 
feedback between the X-ray cooling and the kinetic luminosity of 
the central AGN. 

In Fig. [5] there is also a general trend in the redshift of the 
cluster as the bubble power increases, from bottom left to top right. 
We do not see any high-redshift, low-power bubbles, which is not 
unexpected. At high redshifts these bubbles would be too small to 
detect with the spatial resolution of Chandra. We do, however, see 
some low-redshift clusters with energetic outbursts - Cygnus A and 
Perseus. This is also demonstrated in Fig. [6] There are few large 
bubbles at low redshifts, as a result of the small sampling volume; 
whereas at high redshifts there are few small bubbles as they would 
not clearly be resolved. 

Figs. [5] and [6] imply that we are missing a number of bub- 
bles, mainly at high redshift, that have not as yet been detected. As 
they have not as yet been seen by Chandra, then other methods for 



© 0000 RAS, MNRAS 000, 000-000 



10 Dunn & Fabian 



10 : 





KJ 










• •o 


• 

• 




o 








o 








o 


o o 


A 

O 






o 


o 




o o 




o • 








▲ 


o 






10 43 


10 44 










erg/s) 





Figure 6. The rate of energy loss via X-ray emission from within the cool- 
ing radius of the cluster versus the average bubble size for the cluster. These 
are for the clusters which harbour clear bubbles. The two triangle symbols 
are for M87 and Perseus where L coo i could not be determined from a sin- 
gle ACIS chip. The colour scale is the logarithm of the redshift, blue is low 
redshift, red is high redshift. 



their detection will have to be developed if their effect on the clus- 
ters is going to be analysed in detail, as will have at best 
2 arcsec spatial resolution. As can be seen in Fig. [4] most clusters 
which have short cooling times harbour radio sources, and some 
have clear bubbles. It could be argued that the natural consequence 
of an AGN at the centre of a cluster is that, if it is currently produc- 
ing jets, it produces bubbles and injects energy into the centre of 
the cluster. The GHz radio emission could then be used as a tracer 
of the current size of the bubbles, and from this the energy input of 
the AGN be quantified. 

Another interpretation is that the AGN in high redshift clus- 
ters are undergoing their first and hence explosive outburst, and 
that they will, over time, head towards a steady balance between 
L coo i and Ptmbbie- However, as Cygnus A and Perseus appear in 
very similar locations to the high redshift clusters and Perseus has 
clear evidence for previous outbursts then this explanation is a little 
uncertain. 

The distribution of radii out to which the energy in the bubbles 
can offset the X-ray cooling (see Fig.[7) show that a large number of 
clusters occur around equality with a tail towards "overheating." As 
most clusters are close to rh ca t /'"cool = 1 it is reasonable to assume 
that even out to z ~ 0.4 there is a balance between AGN heating 
and X-ray cooling, though the caveat of missing small bubbles at 
high redshifts should be noted. 



7.1.3 Duty Cycle 

A significant fraction of the clusters in the combined sample have 
a large heating effect on the cluster. There are two clusters in the 
B55 sample which have r^cat /''cool > 2 (A2052 & Cygnus A) and 
there are four in the BCS sample (A115, Hercules A, ZwC12701 
& ZwC13146). As these outbursts are going to arise in clusters in 
which some form of heating is expected, the parent sample of these 
is 34 (20 from B55 and 14 from BCS). Therefore 18 per cent (6/34) 
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Figure 7. The distribution of radii (as a fraction of the cooling radius) out to 
which the energy in the bubbles can offset the X-ray cooling for all clusters 
in the combined sample which have clear bubbles. Hercules A is off the 
right hand side of this figure (r hcat /r coo i = 13.1). 
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Figure 8. The distribution of the ratio of the average to the "maximum" 
size of the bubbles. The uncertainties are estimated using a Monte Carlo 
simulation of the calculations. 



of clusters in this sample have extreme outbursts, and hence the 
duty cycle is also 18 per cent. 



7.1.4 Bubble Sizes 



We repeat the analysis in iDunn & Fabianl feOOfJ) of comparing the 
observed size of the bubbles to the size expected if they offset all the 
X-ray cooling wit hin r coo i. For a detailed description of the method 
see Section 6. 1 of IDunn & FabiarJ J2006). The "maximum" radius 
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of the bubble is defined as 

iW = Mj^M, (4) 

and is compared to the average radius of the bubble i? av = 
\/R\B%,. We assume that an AGN will produce two bubbles, hence 
each will offset L coo i/2. We use Monte Carlo simulations of the 
data to allow estimates to be placed on the uncertainties in the dis- 
tribution^ Only clusters with young bubbles were used for this cal- 
culation. The distribution is shown in Fig.[8]and still shows a peak 
at around i? av /-R max ~ 0.7. 

The peak has no t moved significantly from that found by 
iDunn & Fabianl <2006l) . which is not surprising as only 3 clusters 
(6 bubbles) ha ve been added. The exp lanations for this distribution 
as outlined in IDunn & Fabianl d2006h are therefore all still valid. 
The bubbles with a small value of fl av /-Rmax are small for their 
host cluster, and so are likely to still be growing. Therefore, it is 
natural to state that all the bubbles with R av / i? max > 1 are large, 
i.e. have detached from the core of the cluster and are expanding 
in the reduced pressure of the higher cluster atmosphere. Some of 
these large bubbles may involve AGN i nteractions which a re not 
well described by bubbles, e.g. Hydra A dNulsen et alj |2005). 

As bubbles spend a larger fraction of their lifetime close to 
their full size rather than growing with small radii, it would be ex- 
pected to see the majority of clusters with -R av /-Rm ax ~ 1. To 
shift the peak in the distribution from 0.7 to 1.0 requires 7? max to 
reduce. This could be achieved by reducing L C ooi, but this is an un- 
likely solution. We have already taken a fairly extreme definition 
for the cooling radius, and it would be more sensible to increase 
the cooling radius to move it closer to other definitions of cool core 
clusters. 



7.2 Evolution in the Cluster Samples 

Comparing the number of clusters requiring heating/harbouring ra- 
dio sources between the BCS and B55 samples show that they are 
very similar. The fraction of clusters harbouring radio sources is 
almost exactly the same for each sample. The fraction of clusters 
requiring heating is slightly larger for the BCS sample (33 versus 
25 per cent). This may be a bias towards Chandra observations of 
cooling clusters, but may also be a true increase in the fraction of 
clusters requiring heating. 

By combining the clusters in the BCS and B55 samples we 
have two well defined samples spanning a redshift range of 0.0 ^ 
z ^ 0.37. However the treatment of the two samples has been 
different, and so when investigating any changes in the clusters with 
redshifts we have to redefine the samples so that selection effects 
do no t influence our r esults . 

IDunn & Fabianl d2006h did no analysis of clusters in the B55 
sample which had no bubbles or central radio source, whereas in 
the analysis presented here, all clusters which had Chandra obser- 
vations were studied. To minimise any selection effect we concen- 
trate on clusters with a radio source and/or bubbles. 

We have chosen to concentrate on the temperature, entropy 
and cooling time profiles, and have not split the clusters on whether 
they harbour clear bubbles or radio sources. The profiles of all 71 



clusters are shown in Fig.|9]binned into four redshift bins with equal 
numbers of clusters in each. 

An initial glance at the temperature profiles gives the impres- 
sion that the lowest redshift clusters are cooler, with central temper- 
atures to below 1 keV. However, this is more likely to be a result of 
resolution rather than redshift evolution. The innermost annulus is 
going to correspond to a smaller real radius in the nearby clusters 
as their emission extends over a larger area of the sky and they are 
also going to be brighter. 

The entropy profiles also vary between the lowest and highest 
redshifts, however this is also most likely the result of resolution. 

The cooling time profiles, however, do show some change be- 
tween the redshift bins which is more likely to be the result of evo- 
lution. For a given radius in the cluste^ the lowest redshift clusters 
have a longer cooling time than the highest redshift clusters. How- 
ever, this may still be the result of resolution. In the high redshift 
clusters the inner regions of the cluster may all be merged together 
in a single region. This would reduce the average cooling time of 
this bin. Whereas in the nearby clusters, this region would be split 
into many separate regions with those with shorter cooling times 
occuring at smaller radii. The bulk of the scatter in this plot does 
reduce at higher radii, which is compatible with this argument. 



8 CONCLUSIONS 

We have extended our study of the effect of central AGN in clusters 
to clusters with z > 0.1, using the BCS. There are only six clus- 
ters which harbour clear bubbles, and on average they overcom- 
pensate for the X-ray cooling. These clusters are combined with 
an updated analysis of those in the B55 cluster sample. The aver- 
age distance, as a fraction of the cooling radius, out to which the 
bubbles can offset the X-ray cooling from the combined sample is 
rheat/^coca = 0.92 ± 0.11 for clusters where this distance can be 
reliably determined. Adding in all clusters with bubbles results in 
rheatAcool = 2.10 ±0.40. 

Although it appears as if the AGN bubbles can easily counter- 
act the X-ray cooling within the cooling radius, there is a selection 
effect, especially at the higher redshifts, which means we are miss- 
ing small bubbles in the most distant clusters. 

Using the result from the BCS bubbles, we calculate the ex- 
pected size of bubbles within the clusters which harbour radio 
sources. In three cases (A1068, A1763 and A2204) the predicted 
bubble sizes and observed radio emission are similar. In some of 
the clusters the radio source identified in the NVSS is offset to 
the cluster centre in VLA observations, and so these may be back- 
ground AGN rather than associated with the cluster. A category of 
clusters with cool cores, radio sources and no clear bubbles exists 
in both the B55 and BCS samples. The lack of bubbles is explained 
by low signal-to-noise ratio in the X-ray images of some but not all 
of these objects. 

The distribution of the central cooling time of the clusters 
shows that there are no clusters with bubbles which have t coa \ > 
1.2 Gyr. Also there are only 10/24 clusters with t coo \ > 4.0 Gyr 
which harbour a central radio source, only 5/47 clusters with 
tcooi < 4.0 Gyr have no central radio source. 



8 The values of fi av /fimai were binned 10 4 times using a Gaussian dis- 
tribution for the input uncertainties. The interquartile range in the values for 
the bins is shown by the error bars in Fig. [8] 



9 We note that these profiles have not been scaled to make them truly clus- 
ter independent. They do however align very well without any scaling ap- 
plied 
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Figure 9. The distributions of the combined sample LEFT: Temperature, CENTRE: Entropy and RIGHT: Cooling time. The colour scale shows the redshift of the 
clusters. The redshifts have been binned into four bins with equal numbers of clusters (18) in each bin (blue = z 0.0521, green = 0.0521 < z $C 0.161, 
yellow = 0.161 < z ^ 0.214 and red = z > 0.214). 



We investigated the evolution of cluster parameters with red- 
shift and over the range 0.0 < z < 0.4 we do not find any signifi- 
cant change. The cooling time profiles show some variation but this 
can be explained as being the result of the resolution differences as 
the distances increase. 
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